Uncovering the Hidden Order in URu2Si2 by Impurity Doping 
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We report the use of impurities to probe the hidden order parameter of the strongly correlated 
metal URu2Si2 below the transition temperature To ~ 17.5 K. The nature of this order parameter 
has eluded researchers for more than two decades, but is accompanied by the development of a 
partial gap in the single particle density of states that can be detected through measurements of 
the electronic specific heat and nuclear spin-lattice relaxation rate. We find that impurities in the 
hidden order phase give rise to local patches of antiferromagnetism. An analysis of the coupling 
between the antiferromagnetism and the hidden order reveals that the former is not a competing 
order parameter but rather a parasitic eff'ect of the latter. 

PACS numbers: 76.60.-k, 71.27.+a, 74.62.Dh 



The heavy fermion URu2Si2 has received considerable 
attention because it undergoes a phase transition to a 
state which is poorly understood. The strong interac- 
tions between the U 5f electrons and the delocalized con- 
duction electrons give rise to an enhanced Sommerfeld 
coefficient 7 = 180 mJ/mol K^^ and two phase transi- 
tions at low temperature: the hidden order (HO) transi- 
tion at To 17.5 K gaps approximately 70% of the Fermi 
surface area and a superconducting transition at Tj, ~ 1.4 
K emerges from the remaining charge carriers [l|, |2| . The 
large entropy associated with the HO phase transition 
is suggestive of spin density wave order, yet direct spin 
probes have shown no evidence of intrinsic magnetic or- 
der in pure crystals. Although the HO phase of URu2Si2 
is not itself magnetic, this phase is closely related to an- 
tiferromagnetism (AF) of the U electron spins. Early 
neutron scattering and muon spin rotation (/^SR) studies 
reported a tiny ordered magnetic moment of 0.03 ^jlb /U 
in pure URu2Si2, which led to the concept of small mo- 
ment antiferromag nctism (SMAF) 0,13. However, later 
/iSR and nuclear magnetic resonance (NMR) measure- 
ments tell a quite different story [2,[i]. They reveal an 
inhomogeneous coexistence between small regions of an- 
tiferromagnetic order and hidden order in pure URu2Si2, 
with a relative fraction that tends toward bulk AF under 
pressure 0, H, Substituting Rh for Ru in URu2Si2 
leads to a suppression of the long range hidden order, 
and recent neutron scattering studies revealed large mo- 
ment AF coexisting with the hidden order for large Rh 
concentrations 

In order to investigate the microscopic effects of the 
Rh dopants on the hidden order phase and to characterize 
the emergent AF in U(Rui_2;Rhj.)2Si2 we have measured 
the ^^Si NMR spectrum as a function of temperature and 
Rh concentration. Fig. [1] shows a series of such spectra. 
The resonance frequency of the ^^Si (nuclear spin I = \) 
is given by / = "/Hq{1 + K), where 7 is the gyromagnetic 
ratio of the ^^Si, Hq is the apphed external field (7 T), 
and K is the Knight shift arising from the hyperfine cou- 
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FIG. 1: (color online) ^^Si spectra in U(Rui_iRha;)2Si2 as 
a function of temperature and Rh concentration. The spec- 
tra were obtained by summing the Fourier transforms of Hahn 
echoes for several difi'erent frequencies in a fixed external field 
along the c-direction. (a) Spectra as a function of Rh concen- 
tration 2; at 4 K. Spectra were normalized to the height of cen- 
tral line. The spectra are consistent with commensurate an- 
tiferromagnetic ordering with moments aligned (anti)parallel 
to the c-direction. (b) Temperature dependence of the spec- 
tra for fixed x = 0.02. The antiferromagnetic signature dis- 
appears gradually near 9 K without change of the ordered 
moment. 



pling between the nuclear and electron spins in the solid. 
Aside from a slight suppression of K connected with the 
opening of the gap, there is no visible change of the spec- 
trum at the hidden order transition, Tq. However at a 
lower temperature Tn we find the emergence of two satel- 
lite peaks on either side of the central resonance. These 
satellites arise because of the presence of a static internal 
hyperfine field, i?hf , associated with commensurate anti- 
ferromagnetic order with moments pointing along (001). 
The nuclei resonate in the local field Hq + Hm, where 
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FIG. 2: (color online) The coexistence phase diagram of the 
HO (black square; left axis) and the AF order (red circle; left 
axis) as a function of doping. Blue triangles (right axis) in- 
dicate the magnetic moment of the AF order. The transition 
of the HO order is determined by the specific heat, while the 
induced AF transition is determined by the disappearance of 
the AF peaks. The lines are the calculated To{x) (solid black) 
and fio{x) (dashed blue) using the GL model described in the 
text. The calculated curves are renormalized by the critical 
concentration Xc 0]. INSET: Specific heat over temperature, 
C/T, as a function of temperature T and doping x. There is 
no sign of a bulk phase transition at Tap (a;). 



i^hf = ^MOj Mo is the ordered U spin moment and A is 
the hyperfine coupling. We find that A is unchanged from 
the pure compound (3.6 kOe/zis) Q, and can therefore 
directly measure the antiferromagnetic order parameter, 
M(T,x) ~ /io(T,x), shown in Fig. H 

The spectra in Fig. [T] reveal an inhomogeneous mix- 
ture of antiferromagnetic (satellite peaks) and hidden or- 
der (central peak) regions below T^. We see no specific 
heat anomaly or critical slowing down at T/v, suggesting 
that this transition is not a new thermodynamic phase, 
but rather a crossover to an inhomogeneous coexistence 
Q. The volume fraction of antiferromagnetic domains, 
shown in Fig. [3l varies with both temperature and dop- 
ing. The antiferromagnetic fraction saturates at low tem- 
perature at a maximum of 80% at x = 0.025. These 
observations suggest that the antiferromagnetic domains 
nucleate around the Rh dopants, forming patches with a 
radius ^af on the order of two to three lattice spacings 
at zero temperature. The satellites in the spectrum arise 
from nuclei within these patches of AF, whereas the cen- 
tral resonance arises from nuclei outside. The percolation 
limit is reached at x = Xc where the antiferromagnetic 
patches overlap. 

A priori these results imply that the AF, characterized 
by the order parameter M, competes with the hidden or- 
der, characterized by an order parameter Much like in 
the vortex cores of the cuprates, a competing antiferro- 



magnetic order parameter can emerge in spatial regions 
where the dominant superconducting order parameter is 
suppressed locally [13, EH- Microscopically the impuri- 
ties can create local strains that may stabilize M in the 
vicinity of the Rh However, if this were the case, 

then /io(0, x) should increase with doping and long-range 
antiferromagnetic order would develop above the percola- 
tion threshold at Xc [l3| . For example, a competing order 
parameter is stabilized in Cd doped CeCoIns where anti- 
ferromagnetic droplets are nucleated at Cd dopants and 
long-range order develops when these droplets overlap 
14|. However, in U(Rui_2;Rha:)2Si2 detailed measure- 



ments of the specific heat as a function of both temper- 
ature and doping (see inset of Fig. ^ show no evidence 
of a second phase transition associated with long-range 
antiferromagnetic order, either within the hidden order 
phase or outside the phase when Tq = 0. Furthermore, 
as seen in Figs. [T] and [2] the antiferromagnetic order pa- 
rameter, M(0,a:;) ^ IJ-oi^), vanishes before ^ does. In 
fact, we find that M{0,x) scales with To{x) (Fig. H]), 
suggesting that the AF is controlled by the hidden or- 
der and never exists on its own as true long-range order 
but rather as a parasitic effect within the hidden order 
phase. It is possible that there are in fact two competing 
effects with Rh doping: local strains that stabilize M and 
modifications to the electronic structure from the excess 
carriers introduced by Rh that destabilize both M and \E' 
simultaneously. In this case there is no reason for M{x) 
and ^'(0, a;) to have the same behavior, and the simulta- 
neous disappearance of both order parameters implies an 
unlikely coincidence. 



Recently Elgazzar [15[ has suggested that the HO is a 
dynamic phenomenon in which the Fermi surface is par- 
tially gapped to a commensurate antiferromagnetic state 
that becomes static under pressure. The Rh dopants 
might then serve to pin the local fluctuations of M^, giv- 
ing rise to local static patches. Once again, it is not clear 
why M{x) should track To(x) and the HO is completely 
suppressed when the local patches overlap. It is possible 
that the Rh doping simultaneously pins the fluctuations 
and destabilizes the HO via modifications of the elec- 
tronic structure. As argued above, though, this scenario 
requires an unlikely coincidence. It also is unclear why 
the pinning would take place only within a few lattice 
spacings of the dopant even though the HO is a long 
range phenomenon. Furthermore, we note that there is 
little difference in the temperature dependence of the spin 
lattice relaxation rate measured in the regions of the bulk 
(central peak) versus the AF droplets (satellite peaks), 
suggesting the absence of either a dynamic phenomenon 
or competing order parameter. 

In fact, the observed correlation between the AF or- 
dered moment and the HO gap suggest that the an- 
tiferromagnetic patches are an epiphenomenon that is 
a direct consequence of the local suppression of the 
hidden order in the vicinity of the dopants. We pro- 
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FIG. 3: (color online) The spectral weight of the AF signal 
as a function of temperature and doping. Upper panel: The 
paramagnetic (PM) fraction, defined as the relative intensity 
(area) of the central line. Spectral intensities were corrected 
for the Boltzmann factor, and normalized to the high temper- 
ature values. Lower panel: The fractional area under the AF 
satellites relative to the area of the central PM peak in the 
spectra (Fig. [1]). INSET: The nonmonotonic behavior of the 
AF fraction versus doping at 4 K. Filled squares are the mea- 
sured fraction as shown in lower panel, while empty squares 
are the indirect results of the lost fraction of the PM peak 
in the upper panel. The solid line is the calculated average 
magnetization (scaled AF fraction) discussed in the text. 
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FIG. 4: (color online) {TiT)~^ as a function of temperature 
and doping concentration x. Data are shown for x = (solid 
square, cyan), x = 0.01 (solid diamond, purple), x = 0.015 
(solid down triangle, green), x = 0.02 (solid up triangle, blue), 
X = 0.025 (solid hexagon, red), and x = 0.03 (solid cir- 
cle, black); open symbols correspond to the antiferromagnetic 
satellites at a; = 0.02 (open up triangle, blue), and x = 0.025 
(open hexagon, red). The solid and dotted lines are guides 
to the eye, and the colored arrows indicate To{x). INSET: 
(TiT)"^/^ at T = 4 K versus ATq, revealing the increase in 
A^(0) as the gap is filled by impurity states. 



pose that the antiferromagnetic order is coupled to the 
spatial derivatives of ^'(r). To interpret the results 
we used the Ginzburg-Landau (GL) [1^ free energy 
functional of the combined system that can be writ- 
ten as F[«',M] = Fho + Faf + Fc, with FhoI"^] = 
ai(r-ro)*^ + ^^'i«''* + Ki|V*|2 + y,5(r)1'2^ FAF[yi\ = 
a2|M|2 + ife2|M|4 + ^2 ((VA4)2 + [VMyf + {VA'hf), 
with GL coefficients 01,02,61,62 > 0, and impurity po- 
tential V. The coupling term is Fc[^ , M] = gi^-^jMp + 
52 1 M 1 2 1 V * p + .93 1 M • V * p . The consequences of the first 
coupling term gi have been discussed before 17[, while 



terms 172 and 173 give rise to nucleation of inhomogeneous 
antiferromagnetic order around the impurity site where 
the hidden order is suppressed. Since there is no exper- 
imental evidence for long-range antiferromagnetic order 
in the undoped system at zero pressure, (02,31 > as 
well as 61, 62 > 0), the only way to stabilize a local solu- 
tion of M(r) around an impurity is by demanding that 
92,93 < 0. To simplify our discussion, we consider only 
the coupling term 33 and choose M = (0, 0, M) along the 
applied magnetic field. The eS^ect produced by a 52 term 
would be similar to that of 33 and will be neglected. From 
our analysis it follows that if the hidden order is locally 
suppressed at the Rh dopants, then AF naturally emerges 
in regions near Rh atoms. The length scale for the recov- 
ery of the hidden order, the coherence length £,{T), will 
then determine the spatial extent of the antiferromag- 
netic patches and the percolation threshold then corre- 
sponds to a suppression of the long-range hidden order. 
Assuming for simplicity that the suppression occurs peri- 
odically, then the suppression of the hidden order transi- 
tion temperature will be ATq = Tq—Tq ^ g^Ki/oi, where 
Tq is the transition temperature of the undoped system 
and the wave vector q describes the spatial modulation of 
the hidden order parameter. Thus, to leading order the 
To suppression will be proportional to the square of the 
Rh concentration, as experimentally observed, instead of 
the usual linear dependence for impurity-averaged the- 
ories. Furthermore, as the hidden order 4" is gradually 
suppressed by Rh dopants and Tq diminishes, the induced 
(or parasitic) antiferromagnetic order M will decrease as 
well. If we treat M as as small perturbation to 'f, then 
the maximum value Mq at the impurity site will decrease 
according to Mo(T)2 ^ (2-i|g3|e(T)-2^-o(T)2 - 202)762, 
where the uniform solution of the unperturbed hidden or- 
der is *o(T)^ = -oi(T)/6i. This trend is clearly visible 
in the data in Fig. 2. 

In order to characterize the low energy density of states 
associated with these localized states near the Rh impuri- 
ties, we have measured the nuclear spin lattice relaxation 
rate, Tj~^, as a function of temperature and doping both 
within and outside of the antiferromagnetic patches. As 
seen in Fig. g (TiT)"! - N{0)^ is suppressed below Tq 
because of the development of the partial gap in the den- 
sity of states (DOS) N{0) at the Fermi surface. With in- 
creasing doping, {TiT)~^ increases monotonically within 
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the hidden order phase. This behavior is very similar 
to the effect of impurities in unconventional supercon- 
ductors, suggesting that the Rh impurities induce extra 
states at low energies [isl . [lij . In this case, we expect 
(riT)-i(x) - N^{0,x) ~ (ATo)2(x), consistent with our 
observations (INSET, Fig. S]). Indeed, (TiT)-^ is faster 
at the antiferromagnetic satellites in the spectrum, sug- 
gesting an excess local DOS within the droplets. In the 
case of URu2Si2 there are multiple bands and one expects 
generally two distinct scenarios for the gap to fill up. The 
first corresponds to a gap in the low energy states for all 
bands, in which case the impurity doping would fill up 
the DOS for all the bands. A second possibility is that 
some of the states remain gapless below Tq while others 
develop a full gap, which is consistent with the specific 
heat anomaly. In this case the impurity induces intra- 
gap states in the gap and essentially does not affect the 
DOS of the ungapped states. We expect the latter to 
be realized here. This behavior is also consistent with a 
subsequent superconducting transition observed at lower 
temperatures. To test this scenario one would need to 
observe the DOS in URu2Si2 as a function of Rh in tun- 
neling experiments, like scanning tunneling spectroscopy. 

The fact that the AF is manifest only through the spa- 
tial gradients of the hidden order rules out theories of or- 
bital AF and helicity order 0, HH. Rather, the hidden 
order appears to involve compensated spin polarizations 
on multiple sites such as the triple-spin correlator sce- 
nario or an unconventional multi-band spin density wave 
(2^ . [ist . Similar effects are well known in the study of 
the NMR hyperfine field at the oxygen sites in the doped 
high temperature superconducting cupratcs: Zn or Ni 
impurities substituted at the Cu sites locally perturb the 
staggered antiferromagnetic order of the Cu 3d spins, giv- 
ing rise to finite hyperfine fields at the O sites [2J|. In 
the absence of impurities, the hyperfine field at the O 
site vanishes by symmetry. In U(Rui_^Rh2;)2Si2, gradi- 
ents of the HO parameter may lead to non-cancelation of 
the net spin per U site, giving rise to the static that 
we observe. Our results are consistent with induced mag- 
netism Al being commensurate with the lattice, while the 
hidden order ^' is very likely incommensurate, as was ar- 
gued by Wiebc et al. We point out that this discussion 
implies that Rh doping induces the conversion of HO 
to commensurate AF state within each droplet. If in- 
deed the HO state represents an incommensurate charge 
density wave (CDW), as argued in [i^l then one would 
expect that impurities induce spin dependent scattering 
that converts CDW order into magnetic excitations and 
in addition modifies the momentum of the density wave 
to make it commensurate. Thus it remains a fascinat- 
ing theoretical and experimental puzzle to explain the 
sudden conversion of incommensurate hidden order into 
commensurate antiferromagnetic order in the presence of 



disorder and possibly pressure. 



[1] T. T. M. Palstra, A. A. Menovsky, J. v. d. Berg, A. J. 

Dirkmaat, P. H. Kes, G. J. Nieuwenhuys, and J. A. My- 

dosh, Phys. Rev. Lett. 55, 2727 (1985). 
[2] K. Matsuda, Y. Kofiori, T. Kohara, K. Kuwahara, and 

H. Amitsuka, Phys. Rev. Lett. 87, 087203 (2001). 
[3] C. Brofiolm, J. K. Kjems, W. J. L. Buyers, P. Matthews, 

T. T. M. Palstra, A. A. Menovsky, and J. A. Mydosh, 

Phys. Rev. Lett. 58, 1467 (1987). 
[4] D. E. MacLaughlin, D. W. Cooke, R. H. Heffner, R. L. 

Hutson, M. W. McElfresh, M. E. Schillaci, et al., Phys. 

Rev. B 37, 3153 (1988). 
[5] G. M. Luke, A. Keren, L. P. Le, Y. J. Uemura, W. D. Wu, 

D. Bonn, L. Taillefer, J. D. Garrett, and Y. Omacrnuki, 

Hyperfine Interact. 85, 397 (1994). 
[6] M. Yokoyama, H. Amitsuka, S. Itoh, I. Kawasaki, 

K. Tenya, and H. Yoshizawa, J. Phys. Soc. Jpn. 73, 545 

(2004). 

[7] Y. Kohori, K. Matsuda, and T. Kohara, J. Phys. Soc. 

Jpn. 65, 1083 (1996). 
[8] A. Villaume, F. Bourdarot, E. Hassinger, S. Raymond, 

V. Taufour, D. Aoki, and J. Flouquet, Phys. Rev. B 78, 

012504 (2008). 

[9] In our calculations, we used GL parameters ai = bi — 

b2 = Ki — 1, and 02 = 0.02, K2 = 0.05, with coupling 

coefficients gi = g2 — and 33 = —0.5. 
[10] E. Demler, S. Sachdev, and Y. Zhang, Phys. Rev. Lett. 

87, 067202 (2001). 
[11] S. A. Kivelson, D.-H. Lee, E. Fradkin, and V. Oganesyan, 

Phys. Rev. B 66, 144516 (2002). 
[12] M. Yokoyama, H. Amitsuka, K. Tenya, K. Watanabe, 

S. Kawarazaki, H. Yoshizawa, and J. A. Mydosh, Phys. 

Rev. B 72, 214419 (2005). 
[13] Y. Zhang, E. Demler, and S. Sachdev, Phys. Rev. B 66, 

094501 (2002). 

[14] R. R. Urbano, B.-L. Young, N. J. Curro, J. D. Thompson, 
L. D. Pham, and Z. Fisk, Phys. Rev. Lett. 99, 146402 
(2007). 

[15] S. Elgazzar, J. Rusz, M. Amft, P. M. Oppeneer, and J. A. 

Mydosh, Nat. Mater. 8, 337 (2009), 0809.2887. 
[16] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
[17] N. Shah, P. Chandra, P. Coleman, and J. A. Mydosh, 

Phys. Rev. B 61, 564 (2000). 
[18] A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. 

Phys. 78, 373 (2006). 
[19] S. Ouazi, J. Bobroff, H. AUoul, M. L. Tacon, N. Blan- 

chard, G. Collin, M. H. Julien, M. Horvatic, and 

C. Berthier, Phys. Rev. Lett. 96, 127005 (2006). 
[20] P. Chandra, P. Coleman, J. A. Mydosh, and V. Tripathi, 

Nature 417, 831 (2002). 
[21] C. M. Varma and L. Zhu, Phys. Rev. Lett. 96, 036405 

(2006). 

[22] V. Barzykin and L. P. Gor'kov, Phys. Rev. Lett. 70, 2479 
(1993). 

[23] A. V. Balatsky, A. Chantis, H. P. Dahal, D. Parker, and 
J. X. Zhu, Phys. Rev. B in press (2009), 0903.2570. 

[24] S. Ouazi, J. Bobroff, H. AUoul, and W. A. MacFarlane, 
Phys. Rev. B 70, 104515 (2004). 



